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I. INTRODUCTION 

This research was first proposed in May, 1986, to focus on some of the 
problems encountered in the analysis of postural responses gathered from 
crewmembers, The ultimate driving force behind this line of research was the 
desire to treat, predict, or explain "Space Adaptation Syndrome" (SAS) and 
hence circumvent any adverse effects of space motion sickness on 
crewmember performance. The aim of this project was to develop an easily 
implemented analysis of the transient responses to platform translation that can 
be elicited with a protocol designed to force sensorimotor reorganization, 
utilizing statistically reliable criterion measures. 

This report will present: 1) a summary of the activity that took place in 
each of the three funded years of the project, 2) discussion of experimental 
results and their implications for future research, and 3) a list of presentations 
and publications resulting from this project. Most of the figures and the two 
tables referred to in the text are placed at the end of the body of the report, 
before the appendices. 
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II. SUMMARY OF YEARLY ACTIVITY ON THE GRANT 

Reports were submitted to NASA at the end of the first and the second 
year of the grant period. The original award took place in April, 1988, so the first 
year period covered 4/1/88 - 3/31/89. The first and second year information 
presented here comes directly from the reports submitted: 

First Year 


GRANT ADMINISTRATION 

The starting date for the first year was April 1, 1988. Other dates were 
noteworthy: the NASA postural measurement system was shipped to UMASS 
and arrived August 16. 1988, but the LSI-1 1 system had no terminal for user 
interaction. NASA subsequently shipped a VDT system with keyboard, which 
arrived September 16, 1988. The other event of note was the execution of a 
no-cost extension from April 1, 1989 to June 30, 1989. We received approval 
from our institutional Human Subjects Review Board to perform the experiments 
in August, 1986, and this approval was renewed for 1988-89. 

ACTIVITIES AND ACCOMPLISHMENTS 


Hardware 

The activities during the first year of the grant concentrated on efforts to 
refurbish the LSI-1 1 based NASA -JSC postural measurement system. These 
efforts are chronicled first. 

Upon arrival of the VDT and terminal in 9/88, the system was powered-up 
for the first time. Within fifteen minutes from initial power-up, the dual floppy disk 
(RX01) system spontaneously went into a rapid series of seeks, and 
proceeded to crash. The controller card was visually inspected and re- 
seated, but the floppy drive would not respond. The entire system was 
rendered useless because of the reliance on the disk drives for data 
acquisition and stimulus delivery software. After much digging, two discarded 
LSI-1 1 computers were discovered in moth-balls in the Electrical Engineering 
Department at UMASS. Through December, 1988 and January, 1989 the 
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components from the three LSI-1 1 systems were mixed and matched until a 
final hybrid evolved. The hybrid system included double sided floppy drives 
(RX02), an LSI-1 1/23 CPU, and bypassed the SMU controller from the original 
system (which was suspect). This hybrid system was used to move the posture 
platform and collect data from the digital line scan cameras. Two problems 
still remained: 1) we had no means to inspect the data we acquired because 
the software previously written for the system relied on Tektronix graphics 
capabilities, and our VDT had no graphics capabilities, and 2) we had no 
system or applications software to modify any of the existing software for our 
purposes. These problems were solved by unearthing RT-1 1 v. 5.02 and 
FORTRAN 77 from various sources on campus, in addition to finding emulator 
software for Tektronix graphics terminals to run on IBM compatible PC's. By 
February, 1989, we could finally run the NASA/JSC postural measurement 
system and view the collected data. At that point in time, we feit that with 
some minor alterations to the software to customize our protocols, we were 
ready to start the first phase of data collection. Unfortunately, more serious 
snags developed. 

In the process of running pilot sessions on lab personnel, the 
hardware/software system experienced intermittent failures which took the 
form of either losing control of the moving platform in the midst of a trial, or 
experiencing a hang-up in the software which required system reboot. The 
former occurrence occasionally resulted in high speed platform motions to the 
microswitched cut-off of the servomotor, an event of extreme danger to both 
the subject and the drive screw mechanism. As with the case of software 
hang-up, all previous data was lost and the session had to be restarted. Under 
many environmental conditions, using protocols of varied length (number of 
trials) and waveforms (sinusoidal, transient, and pseudorandom), between the 
two types of failures approximately 50% of the data acquisition was flawed. 

During this time period (1/89-3/89) we were planning the modifications 
needed to the LSI-1 1 data acquisition programs for our needs, The major 
modification was the need to acquire the EMG data at 1000 Hz for four 
channels concurrently. The original software acquired the data at 25 Hz, much 
too slow for our purposes. The limitations of the double-buffering scheme for 
writing data to floppy disk combined with the limited RAM of the LSI-1 1 (128K) 
made the 1000 Hz unattainable. With this realization, we decided to use the LSI- 
1 1 system for stimulus delivery and digital line scan camera data acquisition 
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only, and to acquire the EMG data and platform position data with an IBM 
compatible where the eventual analysis would be performed. We wanted to 
use the digital line scan cameras because they provide a real-time sagittal 
plane body silhouette position; unfortunately, they are mated to the LSI-1 1 bus 
and can not be interfaced directly to the PC without new interface cards. We 
were developing the software on both the LSI-1 1 and the PC to establish this 
stimulus delivery and data acquisition scheme, when the LSI-1 1 refused to boot 
up (6/89). Between the previous intermittent problems and the final difficulty 
booting up the system, I decided to abandon the LSI-11 based approach to 
our experiments at the end of the first year. Luckily, we had been developing 
software for the PC to analyze our data, so that progress was made on that 
front, 

The other hardware related activity scheduled for the first year was the 
specification and purchase of an additional DC torque motor to add an ankle 
rotation degree of freedom to the existing translating platform. The 
calculations of required peak torque indicated that an Industrial Drives TT-2033- 
1010-A motor that produces 100 oz-in. of torque with a maximum speed of 6000 
RPM will be sufficient for our purposes. The motor costs $51 7, and will be 
purchased at the start of the second year. 

Software 

This next section describes the software developments that took place 
throughout the first year, and those that took place at the end of the first year 
when the LSI-1 1 was abandoned. A data analysis package was developed 
to process the EMG signals. This package will accept up to 16 channels of 
analog EMG data (either raw or RMS), perform post-processing if necessary 
(integrate with a given time constant over a given window), and display the 
data (either one channel or multiple channels simultaneously). At this point the 
latency from the onset of platform motion can be manually digitized under 
cursor control or automatically selected and displayed by an algorithm 
adapted from previous work in our lab. The onset and duration of muscle 
activity are determined this way and stored to disk for each EMG channel. This 
software will operate on two different A/D boards: Metrabyte DASH-16 or 
Data Translation DT-2801 . 
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Without the LSI-1 1 system to provide the velocity commands to the 
moving platform, the PC was needed for stimulus generation, The DASH-16 
board has two channels of D/A, so software has been developed to output 
desired velocity commands to the platform. Output from the D/A is 
monopolar, so a bias voltage is used downstream to provide bidirectional 
control. Unfortunately, we discovered that the maximum sample rate to attain 
A/D data concurrently with D/A output under interrupt control was 
approximately 200 Hz, even running QuickBasic 4.0. After an extensive search 
of available products (at a reasonable cost) that would allow us to run 
concurrent A/D and D/A under interrupt control, none were found. Our solution 
to this dilemma is summarized in the section on our finalized data acquisition 
plans. 

When the use of the digital line scan cameras was ruled out, we began 
to develop the software to perform batch analysis of body sway motion 
obtained with a video based system (MotionAnalysis System, Santa Rosa, 

CA). Our system is a 2-D video processor hosted by a SUN 3-160 Workstation, 
and uses two 60 Hz video cameras. Programs have been developed in the 
ExpertVision environment to provide sagittal plane coordinates of the fifth 
metatarsal head, the calcaneus, the lateral malleolus, the knee, the hip, the 
shoulder, the mastoid process, and the elbow. Preliminary experiments 
indicate a resolution of 0.1 cm. 

One advantage of this increase in the amount of body position data 
(from two points with the digital line scan cameras to eight points with the 
video) is that the data will be compatible with a dynamic biomechanical 
model developed over the past year to predict joint reactive moments and 
L5/S1 loading. This will enhance our data analysis and interpretation. 

FINALIZED DATA ACQUISITION PLAN 

The decisions about the final data acquisition system have been made 
and will be implemented in the second year. We will use two 60 Hz video 
cameras and our video-processing system to acquire the body sway data. 
The laboratory IBM-XT will provide the velocity commands to the moving 
platform with the DASH- 16 D/A channel output. The Zenith 386 system 
purchased to analyze the data will have a DT-2801 A/D board installed, and 
will then acquire the platform position and the EMG data at 1200 Hz, The IBM-XT 
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will also output synchronizing signals to both the video processor and the Z-386 
to allow data synchronization. 

Subjects for the first phase of the experiments are being recruited from 
two large undergraduate Exercise Science courses, with the understanding 
that we are looking for a two year availability if possible. 

RFVISFD TIME-LINE FOR SECOND YEAR ACTIVITIES 

Given the delays in establishing a viable data acquisition and stimulus 
delivery system, the second year will be more intensive for actual 
experimentation. Phase I experiments with males will take place in November 
and December, 1989. Phase II experiments with males will be conducted in 
February and March, with the intervening time used to analyze Phase I results 
and to finalize the software for Phase II data acquisition. Meanwhile, the 
females for Phase I experiments will be selected in March, and Phase I 
experiments with females will begin at the end of March (after spring break), in 
April the software for running the Phase III experiments will be completed so 
that Phase III male data can be acquired in May, before the semester ends. 
June of 1990 will be spent analyzing acquired data and preparing the second 
year report. Manuscripts concerning results from Phase I and Phase II will be in 
preparation at this point, 

BUDGET ACTIVITY IN THE FIRST YEAR 

The budget for the first year went according to plan, with some slight 
stretching due to the delay in starting the second year. The equipment budget 
still contains about $2600, but this will be spent early in the second year to 
purchase the second axis DC torque motor (approximately $600) and the A/D 
board and screw terminal for the Zenith 386 (approximately $1500). Many 
supplies purchased during the first year will be carried over into the second 
year (particularly the EMG electrodes and the computer storage media), so 
that the heavy experimental activity of the second year should still remain 
within budget. 



Robert O. Andres, Ph.D. 


8 


SUMMARY 

First year activity concentrated on the development of a reliable stimulus 
delivery and data acquisition system. With the loan of the NASA/JSC postural 
measurement system (providing the recent model moving platform and the 
servoamplifier to add a second degree of freedom) and existing laboratory 
systems (the video-based MotionAnalysis System, the EMG measurement 
systems, and the IBM-XT with the DASH-16 A/D and D/A board), a two degree 
of freedom, dynamic postural measurement system has been assembled for 
essentially the cost of a Zenith 386, a DT-2801 A/D board, and a second DC 
torque motor. Actual hardware and software costs to assemble such a system 
from scratch would easily exceed $200,000; the system now is truly state of the 
art. Second year activity will take advantage of this system to run the 
experiments which we hope to use to design a dynamic posture testing 
protocol to systematically probe the sensorimotor adaptations to postural 
perturbations. 

Year 2 


GRANT ADMINISTRATION 

The starting date for the second year was July 1, 1989. We received 
renewed approval from our institutional Human Subjects Review Board to 
perform the experiments. 

ACTIVITIES AND ACCOMPLISHMENTS 


Hardware 

The data acquisition system described in the First Annual Report was fully 
implemented in the second year. The block diagram in Figure 1 presents the 
final configuration. PC-1 in this figure is the IBM-XT with the Metrabyte DASH-16 
A/D and D/A board. D/A channel 1 drives the platform servo-amplifier, while 
D/A channel 2 sends a trigger signal to both the data acquisition board of PC- 
2 and to each of the VCR's (on audio channel 2, at 400 Hz). PC-2 is the Zenith 
386 with the Data Translations DT-2801 A A/D board. Five channels of data are 
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acquired at 600 Hz for 4 seconds: platform potentiometer, and RMS-EMG's 
from ankle flexors and extensors and knee flexors and extensors. Data are 
saved in packed binary format to conserve disk space. Video records are 
transferred from videotape to digitized marker trajectories by the 
MotionAnalysis System (Santa Rosa, CA) VP-1 10 video processor working with a 
SUN 3-160 based workstation and Expert Vision software. These sagittal plane 
marker trajectories are then downloaded to a PC as ASCII files to allow further 
analysis and to be synchronized with the muscle activity patterns. 



Figure 1 : Block diagram of data acquisition system. 


Software 

A data analysis package was developed to process the EMG signals. 
This package will accept up to 16 channels of analog EMG data (either raw or 
RMS), perform post-processing if necessary (integrate with a given time 
constant over a given window), and display the data (either one channel or 
multiple channels simultaneously), At this point the latency from the onset of 
platform motion can be manually digitized under cursor control or 
automatically selected and displayed by an algorithm adapted from 
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previous work in our lab. The onset and duration of muscle activity are 
determined this way and stored to disk for each EMG channel. This software 
will operate on data collected by two different A/D boards: Metrabyte DASH- 
1 6 or Data Translation DT-2801 . Another analysis approach that has been 
implemented is the calculation of means and standard deviations of RMS- 
EMG's within particular time windows. The first experiment used three time 
windows: before platform motion, during platform motion, and after platform 
motion. 

Experiment 1 

The first experiment was completed in the second year. This experiment 
represented a combination of Phase I testing for both males and females. 
Sixteen subjects (8 male, 8 female) participated in the experiment which 
consisted of three days of testing. Subjects were translated backwards at a 
peak velocity of 25 cm/sec for either 0.3 s or for 1 .5 s. Ten consecutive trials 
were acquired for each duration. Subjects also stood either on the level 
platform or on a slant board tilted backwards by 10 degrees. Therefore each 
session consisted of 40 dynamic trials. The first day was for acclimation 
purposes; the order of tilt and duration presentation for each of the two 
remaining sessions was balanced for order across the subjects. All trials were 
performed with eyes closed, arms folded gently across the chest, and knees 
straight but not locked. 

Analysis of the EMG data is almost complete, but the video data is still 
being processed. Once ASCII files of all the marker trajectories are complete, 
this sway data will be analyzed for settling time (at the shoulder) and will also 
be used as input to a biomechanical model which predicts joint reactive 
forces and moments, particularly at L5/S1. 

RFVISFD TIME-LINE FOR THIRD YEAR ACTIVITIES 

The second experiment will require rearward rotation of the standing 
human, so the first activity will be to add this capability to the platform. The 
servo-motor will be purchased and mechanically coupled to the platform (the 
motor has been selected, as described in the First Annual Report). The 
second experiment will be completed by December, 1990, at which time the 
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results of the first experiment will be in manuscript form. The final 6 months of the 
third year will be spent developing the criteria for an adaptive protocol based 
on Experiments 1 and 2, with Experiment 3 demonstrating the feasibility of the 
approach. At the end of the year specific recommendations for 
implementation of the protocol pre- and post-flight will be made. 

BUDGET ACTIVITY IN THE SECOND YEAR 

The budget for the second year went according to plan, with no major 
deviations from the proposal to report. 

SUMMARY 


Second year activity on the grant concentrated on the actual 
implementation of the reconfigured data acquisition system, culminating in the 
completion of the first major experiment aimed at developing an adaptive 
protocol for astronaut screening. Third yedr activities will complete the stimulus 
delivery system by adding rotation to the platform, and then complete the 
experiments needed to recommend a new protocol which will assess the 
standing human's capability to reorganize their sensorimotor behavior. 

Year 3 


GRANT ADMINISTRATION 

The activity summarized here under Year 3 actually covers the final year 
of grant funding and a one year no-cost extension on the project. One item of 
note was that not all of the money was spent from the budget; money was left 
in to support the construction of the rotating addition to the translating platform 
because of a lack of access to facilities that could machine metal to the 
required tolerances within our department. To have the work done in another 
department, an excessive delay was unavoidable because of staffing 
cutbacks at UMASS and departmental priorities. Hence our efforts were 
concentrated on the analysis and interpretation of the data we had acquired 
with the translating platform. 
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Figure 2: CV's for a 0.3 s translation with no tilt 



Gastrocnemius Tibialis Anterior Hamstrings Quadriceps 


Figure 3: CV's for a 0.3 s translation with rearward 

tilt 



Gastrocnemius Tibialis Anterior Hamstrings Quadriceps 
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Figure 6: Average RMS-EMG's for a 0.3 s translation with no 

tilt 
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Figure 7: Average RMS-EMG's for a 0.3 s translation with 

rearward tilt 



Gastrocnemius Tibialis Hamstrings Quadriceps 

Anterior 
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Figure 8: Average RMS-EMG's for a 1.5 s translation with no 

tilt 
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Figure 9: Average RMS-EMG's for a 1.5 s translation with 

rearward tilt 
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Figure 10: Pre-translation CV's 
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Figure 1 1 : Per-translation CV’s 



Gastrocnemius Tibialis Hamstrings Quadriceps 

Anterior 
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Figure 12: Post-translation CV's 



Gastrocnemius Tibialis Hamstrings Quadriceps 

Anterior 
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Figure 13: Pre-translation Average RMS-EMG's 
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Figure 14: Per-translation Average RMS-EMG's 



■ 0.3 s, no tilt 
□ 0.3 s, tilted 
H 1.5 s, no tilt 
S 1.5 s, tilted 
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Figure 15: Post-translation Average RMS-EMG's 
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ED 0.3 s, tilted 
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Window 

Variable 

Greater on day: 

Pre-translation 

MEANG 

3 


CVTA 

2 


MEANTA 

3 


MEANH 

3 




Per-translation 

MEANG 

3 


CVTA 

2 


MEANTA 

3 


CVH 

2 


MEANH 

3 




Post-translation 

CVG 

2 


MEANG 

3 


CVTA 

2 


MEANTA 

3 


CVH 

2 


CVQ 

3 


Table 1 : ANOVA analysis results with post hoc tests for significant day effects. 
CV=coefficient of variation, MEAN=average RMS-EMG value scaled to 
the maximum for that muscle on that day, G=gastrocnemius, TA=tibialis 
anterior, H=hamstrlngs, Q=quadrlceps. 
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Window 

Variable 

Effects 





Gender 

Duration 

Tilt 






Pre-translation 

CVG 

M>F 


Y>N 


MEANG 



N>Y 


CVTA 

F>M 

L>S 

Y>N 


MEANTA 

F>M 

L>S 



CVH 

F>M 

L>S 



MEAN H 



N>Y 


CVQ 

F>M 

L>S 



MEAN Q 



N>Y 






Per-translation 

CVG 

M>F 


Y>N 


MEANG 

F>M 

S>L 

N>Y 


CVTA 

F>M 

S>L 

Y>N 


MEANTA 

F>M 

L>S 

Y>N 


CVH 

M>F 

S>L 

Y>N 


MEANH 

M>F 

L>S 

Y>N 


CVQ 

M>F 

S>L 

Y>N 


MEAN Q 

F>M 

S>L 

N>Y 






Post-translation 

CVG 

M>F 

L>S 

Y>N 


MEANG 

M>F 


N>Y 


CVTA 


L>S 

Y>N 


MEANTA 

F>M 


Y>N 


CVH 


L>S 

Y>N 


MEANH 

F>M 

L>S 



CVQ 


L>S 



MEAN Q 

M>F 

S>L 



Table 2: ANOVA analysis results with post hoc tests to indicate significant 

gender, tilt, and stimulus duration effects. CV=coefficient of variation, 
MEAN=average RMS-EMG value scaled to the maximum for that muscle 
on that day, G=gastrocnemius, TA=tibialis anterior, H=hamstrings, 
Q=quadriceps, M=male, F=female, L=long duration stimulus (1.5 s), 

S=short duration stimulus (0.3 s), Y=tilted platform (10 degrees of 
dorsiflexion, and N=no tilt. 
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APPENDIX A: 

Manuscript from the presentation at the Xlth International Symposiom of the 
Society for Postural and Gait Research, Portland, OR, May 24-27, 1992. 



Posture and Gait: 

Control Mechanisms, 1992 

Xlth International Symposium of the Society for 
Postural and Gait Research 
Portland, May 24-27, 1992 
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The subjects were required to track a continuously moving target dis- 
played on a monitor while standing on a force plate. The monitor was setted 
up at 1. 5m ahead of subject and the presented target spot was moved recipro- 
cally and slowly on vertical line and CFP signal spot. The velocity of the 
target movement was decided as 0.05, 0.10, and 0.15Hz with a triangular wave 
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APPENDIX B: 


Table B-1 : The calculation of intraclass correlation coefficients. 



Intraclass.nasa 2/5/92 




28 C4V 1 _ 2 2 0.025383 0.0113613 0.026483 0.00140217 0.000055 0.00070109 0.00056807 

29 ~ C5V 1 2 2 0,00032 0.002262 0,1442447 -0.0001942 0.00719624 -0,0000971 0.0001131 

30 ~ AV2G ~ 1 2 2 776.79 27.17 500.18 74.962 -13.8305 37.481 1.3585 

31 AV3G 1 2 2 278.22 32.445 232.03 24.5775 -2.3095 12.28875 1 .62225 
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APPENDIX C: 


Manuscript from the Journal of Biomechanics that presents the dynamic 
biomechanical model that has been modified for this project. 
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VALIDATION OF A BIODYNAMIC MODEL OF PUSHING 

AND PULLING 

Robert O. Andres* and Don B. Chaffin t 

* Department of Exercise Science, 21 Totman Building, University of Massachusetts, Amherst, 

MA 01003, U.S.A. and f Center for Ergonomics, The University of Michigan, Ann Arbor, Michigan, 

U.S.A. 

Abstract — Pushing and pulling during manual material handling can increase the compressive forces on the 
lumbar disc region while creating high shear forces at the shoe-floor interface. A sagittal plane dynamic 
model derived from previous biomechanical models was developed to predict L5/S1 compressive force and 
required coefficients of friction during dynamic cart pushing and pulling. Before these predictions could be 
interpreted, however, it was necessary to validate model predictions against independently measured values 
of comparable quantities. This experiment used subjects of disparate stature and body mass, while task 
factors such as cart resistance and walking speed were varied. Predicted ground reaction forces were 
compared with those measured by a force platform, with correlations up to 0.67. Predicted erector spinae 
and rectus abdominus muscle forces were compared with muscle forces derived from RMS-EMGs of the 
respective muscle groups, using a static force build-up regression relationship to transform the dynamic 
RMS-EMGs to trunk muscle forces. Although correlations were low, this was attributed in part to the use of 
surface EMG on subjects of widely varied body mass. The biodynamic model holds promise as a tool for 
analysis of actual industrial pushing and pulling tasks, when carefully applied. 


INTRODUCTION 

Carts of various sizes, weights and configurations are 
frequently pushed or pulled manually in many in- 
dustrial situations. Among these are the tyre manufac- 
turing industry, the fiberglass manufacturing industry, 
commercial laundries, and the airline industries. A 
large proportion of these tasks involve the worker 
pulling on doors and hoses, pushing on carts, and 
generally attempting a task which imparts a high 
shearing force to the feet, with resulting slip and fall 
injuries (Safety Sciences, 1977). It is estimated that 
over 20% of the worker’s compensation costs each 
year are a result of fall or slip related injuries in the 
U.S. (Szymusiac and Ryan, 1982). A study concerning 
a large manufacturing operation in England reported 
that 36-45% of back pain was caused by a slip or fall 
(Manning, 1983). These statistics paint a dangerous 
picture for workers involved with pushing or pulling 
tasks, which increase the risks of slipping and falling or 
overexerting the back. 

Several investigators have addressed the per- 
formance aspects of push/pull tasks. Dempster (1958) 
studied static pull forces during standing. Kroemer 
measured maximal isometric pushing forces in 65 
different positions and assessed the effects of varied 
foot friction during pushing (Kroemer, 1969, 1971). 
Ayoub and McDaniel (1974) studied the loading of the 
lumbar spine during static pushes and pulls against a 
wall, in various body postures. One Swedish study 
(Winkel, 1983) concerned the manual handling of food 
and beverage carts on wide-body airplanes and meas- 
ured only the hand forces exerted on stationary carts. 
Several recommendations were made about cart con- 
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figuration and loading as a result of this study. How- 
ever, none of these studies considered the dynamic 
case where the worker moves during the push or pull 
task. 

Pushing and pulling hand forces have been meas- 
ured while the subject walked on a treadmill (Snook ex 
ai, 1970) with different handle heights and adjustable 
treadmill resistance. Strindberg and Peterson (1972) 
used psychophysical methods to study force percep- 
tion while pushing trolleys. These studies began to 
approach more realistic dynamic simulations of actual 
industrial situations. A German group studied the 
load on the spine during the transport of dustbins 
(Jager et al. y 1984). These authors utilized a simple 
static model of L5/S1 torques, and they also measured 
the EMG activity of back, leg, and hand muscles. 
However, no validation of their model was offered, 
and the EMG information for the back muscles was 
not compared with the L5/S1 torque predictions. 

The doctoral research of Lee (Lee, 1982; Lee et ai, 
1989) formulated a dynamic biomechanical model of 
cart pushing and pulling. The. inputs to the model 
included subject anthropometry, body postures dur- 
ing dynamic tasks, and hand forces exerted on the cart 
handle. The specific model predictions were horizon- 
tal and vertical foot forces and gross torso muscle and 
vertebral column loadings when pushing or pulling. 
Laboratory validation of the model took place with a 
cart simulator. Dynamic foot forces, hand forces re- 
quired to move the cart, body motions at various 
speeds, and back musde actions were measured while 
six subjects pushed or pulled the cart simulator. Foot 
force predictions were compared to measured (by 
force platform) foot forces, while predicted torso 
muscle forces were converted to ‘equivalent’ integ- 
rated electromyograms (IEMGs) and then compared 
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with the measured IEMGs. These validation experi- 
ments were performed at three handle heights (66, 109, 
152 cm) and three cart resistances (light, medium, 
heavy). Two speeds (approximately 1.8 and 
3.6 km h" 1 ) were used at only the middle handle 
height. The four male and two female subjects ranged 
in stature from 161 to 175 cm, and the body mass 
range was 50-80 kg. 

Lee et al. (1989) reported the results of static calib- 
rations performed for the dynamic experiment, but the 
dynamic results have not been published. The purpose 
of our investigation was to validate a model (based on 
Lee’s model) in the laboratory with a broader anthro- 
pometric population. 


METHODS 

Model development 

The biomechanical model used for this study was 
based on that used by Lee (1982). The mathematical 
relationships and assumptions which allow the predic- 
tion of the desired outputs (foot reactive forces and 
L5/S1 compressive forces) will be described. There 
were several intermediate steps before the predictions 
were made; these will also be examined. 

Two different sets of link lengths were calculated. 
The distance between the joint markers (LEDs), in the 
sagittal plane, was measured by the position detection 
system. The Y and Z coordinates were calibrated first, 
corrected for distortion, and smoothed with a digital 
filter with optimal cut-off frequencies [determined for 
each data record as described by Jackson (1979)] 
before calculating relative link lengths (link lengths in 
the sagittal plane, possibly foreshortened if the link 
was out of plane). Absolute link lengths were based on 
ratios to total stature (Drillis and Contini, 1966). 

The position of L5/S1 was calculated because there 
was no marker on that location. These calculations 
were based on Chaffin and Andersson ( 1984, Chap. 6). 
The link center of gravity (CG) positions were calcu- 
lated with respect to relative link length. Ratios based 
on Dempster (1955) were used. 

Link masses were calculated as a ratio to total body 
mass. A major assumption made at this point was that 
the left arm position was the same as the right arm, as 
viewed from the sagittal plane. When markers were 
placed on the inside of the left arm, they were never 
detected by the position detection system because they 
were blocked by the right arm. Therefore, the masses 
of the forearm and hand link and the upper limb link 
were doubled to represent both arms acting in the 
same sagittal plane position. The ratios were again 
from Dempster (1955). 

Joint angles were calculated from joint marker 
coordinates with respect to the horizontal. The joint 
angles correspond to the angle from the horizontal of 
respective links (i.e. joint angle No. 1 was the angle 
from the horizontal of link No. 1, the forearm and 
hand; see Fig. 1). 



LF, RF Z 



Fig. 1. Free body diagram of double support stance with 
joint angles used by the model and locations of joint center 
markers. Variables and resulting equilibrium equations are 
given in the text. 


The whole body CG position was calculated using 
the segmental moments method (Miller and Nelson, 
1976). The location of the upper body CG was calcu- 
lated with the same approach using only links above 
L5/S1 (hands, forearms, upper arms, and trunk from 
shoulder to L5/SI). 

Mass moments of inertia were calculated for each 
link as the product of the link mass and K 2 y where K 
was the product of relative link length and a constant 
representing the radius of gyration (Chandler et al ., 
1975; Plagenhoef, 1966). 

Link CG, whole body CG, upper body CG, and 
joint angle data were differentiated twice so that the 
respective accelerations could be used for calculation 
of inertial forces. The derivatives were calculated with 
a finite-impulse-response digital recursive filter (Op- 
penheim and Schafer, 1975; Lanshammar, 1982). Opti- 
mal filter coefficients were derived for these data and 
subsequently applied uniformly. 

The linear and rotational inertial forces resist body 
linear and angular accelerations. These inertial forces 
were calculated for each link as the negative of the 
product of the link mass and the linear acceleration in 
Y and Z directions. Whole body inertial forces acting 
at the whole body CG were calculated as the negative 
of the product of the whole body mass and the Y and Z 
accelerations of the whole body CG. Since there was 
no overall body angular acceleration, the whole body 
rotational inertial force was calculated as the sum of 
the individual link rotational inertial forces. 
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Calculation of ground reaction forces at the feet 

Two different situations exist during normal walk- 
ing gait: a single support phase, when only one foot 
contacts the support surface while the other foot 
swings through to its next placement; and a double 
support phase, when both feet are in contact with the 
support surface (see Fig. 1 ). Foot force calculations are 
detailed below for each case. 

Single support . Assume: 

(1) two arms act as one (sum Y and Z forces from 
the separate handles); 

(2) quasi-static equilibrium, 

R F y = H y + B y> 

ZF x =0 = RF x -H x -BW-B x , 

RF t = H z + BW+B t , 

where R F y is the Y reactive force of the right foot; R F x 
is the Z reactive force of the right foot; H, is the total Y 
hand force; //, is the total Z hand force; B W is the 
total body weight; B v is the )' hotly inertial force; und 
B x is the Z body inertial force. 

Moments at the foot during single support were not 
used by the model and were not needed because there 
were only two unknowns with two equations. 

Double support . Assume: 

(1) quasi-static equilibrium; 

(2) two arms act as one; 

(3) moment arms from heel marker for pushing, 
from toe nwker of rear fool for pulling, 

ZF y = 0 = LF y + RF y -H y -B y , 

ZF t = 0 = LF x + RF x -H x -B x -BW. 

Left foot back . 

IM 1 = 0= —RF x *(DF y )+(B x + BW)* DCG, 
-B y *(DCG s )-H y + (DH t ) 

+ //,*(D//,) + B r 

= DCG y * (B x + B W) + DH y * (H t ) + B r 
— DF y + (R F z ) — DCG t + (B y ) 

— DH t + (H y \ 

RF X = (DCG y * (B x + BW) + DH y * (H x ) + B t 
- DCG t * (B y ) - DH X * (H y ))/D F r 
LF x = H t + B x + BW-RF„ 

where L F y is the Y reactive force of the left foot; L F x is 
the Z reactive force of the left foot; DCG y is the Y 
distance from the rear heel to the whole body center of 
gravity; D F y is the Y distance from the rear heel to the 
Z foot force of other foot; DH y is the Y distance from 
the rear heel to the handle; DCG X is the Z distance 
from the floor to the whole body center of gravity; DH X 
is the Z distance from the floor to the handle; and B r is 
the rotational body inertial force. 

Similar equations result when the right foot is back 
for the other portion of double support during right 


foot stance, with appropriate exchange of right and left 
foot reactive forces. At this stage, the system is indeter- 
minate because we have three equations ( M „ F r and 
F x ) and four unknowns ( LF X , RF„ LF r and R F y ). The 
solution requires another equation to become deter- 
minate; the model assumed that the friction utilization 
under both feet was the same, and hence the ratios of 
horizontal (y) to vertical (r) foot forces at each foot 
were equal (Lee, 1982). 


Calculation of trunk muscle forces 

The calculation of L5/S1 compressive forces and the 
muscle forces contributing to these forces begins with 
the calculation of abdominal pressure, because this 
pressure counteracts some of the contraction force of 
the erector spinae muscles (see Fig. 2). An empirical 
prediction of abdominal pressure was performed using 
previously reported equations (Lee et ai , 1989; Chaffin 
and Andersson, 1984) derived from work done by 
Morris et at. (1961). 

The moment arm lit which F Anh nets hus been 
assumed by Chaffin (1975) to vary as the sine of hip 
angle, with an erect position having a moment arm of 
7 cm, increasing to about 1 5 cm when stooped over at 
0 h = 9O° from vertical (where = the angle from the 
hip-to-shoulder link to vertical). The argument that 
F ABD acts parallel to the compressive force on L5/S1 
was presented by Chaffin and Andersson (1984). The 
line of action of rectus abdominus has also been 
parallel to the compressive force on L5/S1 in other 
studies (Schultz und Andersson, 1981; Chullln and 
Andersson, 1984). This model assumes that all muscle 
forces act normal to the shear force to create com- 
pression only. Reactive shearing forces are then pro- 
duced by lumbar facet joints, as described in Chaffin 
and Andersson (1984). 

The following equations were used by the model to 
calculate back and muscle forces (see Fig. 2): 

F c = ESMF + RAMF— F ABD 
+ sin(a) + (BW m + H x + U B x ) 

+ cos (a) * ( U B y + H f \ 

F s = cos (a) * (B W u + H x + V B x ) + sin (a) * ( U B y + H y \ 

where F c is the L5/S1 compressive force; F s is the 
L5/S1 shear force due to external forces only (assum- 
ing that all muscles act to create compression only); 

ESMF is the erector spinae muscle force (when res- 
ultant moment at L5/S1 was negative), calculated as 
the resultant moment at L5/S1 divided by the moment 
arm [0.06 m for males and females, based on Kumar 
(1988)]; RAMF is the rectus abdominus muscle force 
(when resultant L5/S1 moment was positive), calcu- 
lated as the resultant moment divided by the moment 
arm [0.10 m for males and females, from Kumar 
(1988)]; F ABD is the abdominal force due to intra- 
abdominal pressure; a is the angle from horizontal to 
L5/S1 -shoulder link; BW U is the body weight above 
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Fig. 2. Free body diagram of the body above the hips, used 
to calculate L5/S1 shear and compressive forces and erector 
spinae and rectus abdominus muscle forces. All trunk forces 
are acting parallel or normal to the compressive force. (Note: 
it is assumed these muscles act to provide only compression 
forces on the spinal discs.) 


L5/S1; UB X is the upper body Z inertial force; and U B y 
is the upper body Y inertial force. 

Static calibration experiments 

Linear regressions were formed between the predic- 
ted torso muscle forces (erector spinae muscle forces 
[ESMF s] and rectus abdominus muscle forces 
[ RAMFs ]) and the measured RMS-EMGs from both 
the right and the left groups of these muscles during a 
gradual increase of push or pull forces. The subject 
either pushed or pulled in a free static posture. Starting 
with no exerted hand force, the subject built up to a 
near maximum push or pull force in a 5 s period. The 
regression coefficients were used subsequently to 
transform dynamic RMS-EMG values to muscle force 


values to provide a basis for comparison to predicted 
dynamic muscle force values. 

Experimental design and statistical analysis 

A fractional factorial design was used to examine 
the accuracy and precision of the model predictions as 
a function of the independent variables, described 
below. In particular, 20 subjects were selected accord- 
ing to a design which emphasized anthropometric 
extremes ( < 20th percentile, > 80th percentile for both 
stature and body weight) and included subjects of 50th 
percentile stature and weight, with two males and two 
females in each of the five cells (see Table 1). Healthy 
young subjects (18-31 yr) volunteered for this experi- 
ment. After signing an informed consent form, con- 
sistent with University policy, they were weighed and 
their heights were measured. 

Model validation involved three types of analyses: 

(i) correlations of measured vs predicted variables; 

(ii) ANOVA analysis of measured and predicted vari- 
ables; and (iii) ANOVA analysis of residuals formed by 
taking the difference between measured and predicted 
values. SchefTe multiple comparison tests were per- 
formed post hoc , when appropriate. 

Linear statistical models were formed for the 
ANOVA analysis with the following factors (and 
corresponding levels): subject (a = 20), sex (male, fe- 
male), height (short, average, tall), weight (light, aver- 
age, heavy), cart resistance (low, high), walking speed 
(60, 100 steps min -1 ), and direction of resistance 
(push, pull). These factors will be abbreviated SN, SX, 
HT, WT, CR, WS, and DR, respectively. Only the 
main effects are presented in this report; first order 
interactions were not significant. 


Table 1. Subject stature and weight descriptions for experimental design to emphasize 
anthropometric extremes 




Weight (percentile) 




<20 

50 

>80 

Males 

Height (percentile) 

<20 

0.170 m 
626.0 N 

— 

1.75 m 
885.5 N 


50 

I 

1.78 m 
712.0 N 

— 


>80 

1.85 m 
689.0 N 

— 

1.84 m 
878.0 N 

Females 

Height (percentile) 

<20 

1.54 m 
438.0 N 

— 

1.58 m 
845.0 N 


50 

— 

1.62 m 
572.5 N 

— 


>80 

1.72 m 
569.5 N 

— 

1.76 m 
689.0 N 


Stature and weight percentiles were derived from National Health Survey (1965). There 
were two subjects in each cell, so only the mean values are given. 
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Analytical methods 

The differences between independently measured or 
derived variables and the same variables predicted by 
the biodynamic model provided the most critical 
assessment of model performance. Beyond visual com- 
parisons, one technique to assess differences between 
predicted and measured quantities was to compare the 
mean values of the predicted and measured dependent 
variables for each of the two time-windows (single and 
total support period for the right leg). 

Another method of comparing measured and pre- 
dicted quantities involved the formation of residuals. 
The GRF residuals were calculated on a point by 
point basis throughout the support phase by summing 
the model predicted force exerted on the foot by the 
ground and the measured force exerted on the ground 
by the foot. This was equivalent to subtracting meas- 
ured values of force exerted on the foot by the ground 
from the model predicted forces. These residuals were 
actually intra-subject comparisons. Interpretation of 
these results requires some explanation: the average 


represents the sum of the differences between predic- 
ted and measured values divided by the number of 
differences taken. These residual parameters were 
submitted subsequently to ANOVA analysis to deter- 
mine which subject and task factors contributed signi- 
ficantly to the errors in the predictions. 


Data acquisition hardware 
The cart simulator had handles 0.5 m apart and 
oriented so that they were horizontal; hence the 
subject gripped each handle with the hands prone. The 
simulator travelled on teflon bushings over aluminum 
rails, while the cart resistance came from a strap which 
passed over a variable number of dowels which were 
affixed to the bottom of the simulator. By changing the 
tension on this strap, cart resistance varied from 88 to 
128 N (the horizontal force necessary to keep the cart 
moving at approximately 0.5 ms l ). The electronics 
for the portable handles and the EMGs were carried 
by the cart (see Fig. 3). 
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Approximated sagittal plane joint center coordina- 
tes were sensed by a single camera SELSPOT system 
(Selcom Selective Electronic Inc., Valdese, North Car- 
olina, U.S.A.) sampling at 50 Hz. The 10 LEDs were 
placed on the subject as shown in Fig. 1. To prevent 
distortion the following conditions were placed on 
LED position detection: (1) the data window was in 
the central portion of the camera viewing field; (2) lens 
distortion was compensated for by calibration; and (3) 
reflection problems were minimized by keeping the 
reflective laboratory floor out of the viewing window. 
Hand forces were measured with portable handles 
which used columns instrumented with strain gauges 
to detect the horizontal and vertical components of 
force. The force platform was a Kistler force platform 
(Kistler Instrumente AG, Winterthur, Switzerland, 
model 9231 A) with six Kistler charge amplifiers (5001) 
and a Kistler central control unit (5671 A). The surface 
of the force platform was covered with painted ply- 
wood, as was the rest of the runway. The static 
coefficient of friction (CO F) of the surface when rubber 
shoes were worn was approximately 0.7, while the 
dynamic COF was about 0.6 (see Andres et at., 1984). 

The myoelectrical activity of the torso muscles was 
detected by bipolar surface electrodes (Hewlett Pack- 
ard, Andover, Massachusetts, U.S.A., model 14445C 
disposable electrodes). The electrodes were placed 
3 cm lateral to the midline on both sides of the spine at 
the L2 and L3 level, about 5 cm apart so that the 
activity of one side vs the other could be observed. The 
rectus abdominus activity was recorded by placing the 


implies that predictions during double support (that 
portion of the total right leg support period excluding 
single support) were not as valid as those during single 
support. 

Another analysis of the relationships depicts the 
grand mean of the average value of the variable during 
single support by subject body weight category for 
each direction of exertion. The greatest discrepancy 
emerged for the horizontal GRFs during pushing, 
with consistent model overprediction. Model over- 
prediction was also apparent for the vertical GRFs 
during pushing. 

Another parameter selected for analysis was the 
maximum value of the variable within the time- 
window, because the direct comparison between 
measured and predicted values possible with the 
GRFs may have been more sensitive to changes in 
maximum values than to average values (see Table 2). 
Notice again that the slope of the relationship was 
closest to unity (0.7) for the vertical GRF during single 
support. Less of the variability in the data was ex- 
plained by these regressions (using maximum values) 
compared with the average values shown earlier in 
Table 2. 

Comparisons of the trunk muscle forces derived 
from RMS-EMGs with model predicted trunk muscle 
forces indicated that even less of the variability in the 
data was explained by the simple linear regression 
model. Most of these regressions were not significant 
at the 0.05 level. The following section will present 


electrodes in a similar manner on the abdomen (cen- 
tered 3 cm lateral to the linea alba, above the navel). 
The electrode signals were sent to preamplifiers in 
a small box attached to the subject’s belt (input 
impedance = 10* £1, common mode rejection ratio 
>120 db), then to a custom amplifier which con- 
verted the raw EMGs to derived RMS values (with a 
time constant of 55 ms) for the four channels of infor- 
mation (left and right erector spinae, and left and 
right rectus abdominus). 

RESULTS 

Figure 4 displays an example of measured vs model 
predicted results. Quantitative results follow from 
correlation, ANOVA, and residual analyses. 

Correlation analysis of measured vs model predicted 
variables 

Both foot force and trunk force variables were 
compared by forming a simple linear regression be- 
tween measured and predicted values. Comparisons 
were made for the average values of the variables over 


Tabic 2. Average and maximum GRF correlation results 
with linear regression coefficients 

Variable R-square * Intercept Slope 

Average GRFs 


Single support 

RF V 

0.63 

— 10.9 


y 

RF S 

0.66 

27 

0.93 

COF 

0.1 

0.15 

0.1 

Total support 

RF, 

0.67 

-12.9 

0.55 

y 

RF X 

0.12 

357 

0.26 

COF 

NS 

— 

— 


Maximum GRFs 


Single support 

R F v 

0.45 

14.3 

U.J 

y 

RF Z 

0.45 

164.4 

0.78 

COF 

0.06 

0.27 

0.05 

Total support 

RF ¥ 

0.32 

55.4 

0.18 

y 

RF t 

0.21 

435 

0.42 

COF 

NS 

— 

— 


both single and total right leg support phases of the 
gait cycle (see Table 2). Notice that the vertical foot 
force during single support had the slope parameter 
closest to unity. However, for total right leg support 
the predicted average vertical foot force did not correl- 
ate nearly as well with its measured counterpart. This 


RF y ~ horizontal GRF, R F x = vertical GRF, CO F = ratio 
of RF y to RF X . Notice that the correlations during single 
support exceeded those from the total support phase, demon- 
strating the improved performance of the model during single 
support. Also note that the slope parameters were closest to 
unity during single support with the vertical GRFs. 
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Horizontal GRF, pulling 



Time (sec) 


more detailed comparisons of derived and predicted 
trunk muscle forces by performing ANOVA analysis 
with subject and task factors included in the statistical 
models. 

ANOVA analysis of model predicted variables 

The biodynamic model predicted variables and the 
measured variables were analyzed with the inclusion 
of subject (SN, SX, HT, WT) and task (CR, WS, DR) 
factors. Predicted values exceeded measured values for 
both horizontal and vertical GRFs, whether for single 
or total right leg support (as seen by comparing grand 
means in Table 3). Horizontal GRFs were dependent 


always on direction of exertion, while all vertical 
GRFs were affected by individual subject factors. Only 
the mean measured and predicted vertical GRFs 
during the total right leg support phase had the same 
factors in the ANOVA analysis. 

Predicted trunk muscle forces were subjected to the 
same analysis (see Table 3). The model underpredicted 
the ESMFs , but the predicted RAMFs fell between the 
derived values for the right and left side. Since a wide 
range of subject anthropometries were purposefully 
used in these experiments, similar analyses were per- 
formed separately for each of the three weight cat- 
egories (see Table 4). One noticeable difference be- 
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Tabic 3. Average foot and trunk muscle force ANOVA summary 
Variable R-square Mean Subject factors Task factors 


Foot forces 


Single support 

Meas. RF f 

0.84 

-12.18 

— 

DR 

Pred. RF y 

0.69 

1.49 

SN 

DR 

Meas. R F t 

0.66 

582.15 

SN, SX, HT 

DR 

Pred. R F, 

0.86 

609.38 

SN 

— 

Meas. COF 

0.14 

1.15 

— 

— 

Pred. COF 

0.14 

0.29 

— 


Total support 

Meas. R F y 

0.88 

-12.82 

HT 

WS, DR 

Pred. R F, 

0.69 

1.93 

SN 

DR 

Meas. R F t 

0.8 

490.71 

SN, SX, HT 

— 

Pred. RF t 

0.7 

518.93 

SN, SX, HT 

— 

Meas. CO F 

0.18 

0.76 

— 

— 

Pred. COF 

0.29 

0.22 

Trunk forces 

CR, WS 


Single support 

ESMF 

0.39 

561.62 

SN 

CR 

FML 

0.73 

765.2 

SN 

DR 

FMR 

0.71 

858.06 

SN 

DR 

RAMF 

0.32 

71.01 

SN 


FRL 

0.36 

35.48 

SN 


FRR 

0.51 

88.59 

SN 

DR 

Total support 

ESMF 

0.4 

511.23 

SN 

CR 

FML 

0.74 

742.29 

SN 

DR 

FMR 

0.69 

860.35 

SN 

DR 

RAMF 

0.28 

85.79 

SN 

DR 

FRL 

0.36 

32.94 

SN 


FRR 

0.51 

88.06 

SN 

DR 


‘Pred.’ before the variable denotes the predicted value, and *Meas.’ denotes the measured value, RF y 
=* horizontal GRF, RF S « vertical GRF, CO F = ratio of R F y to R F r ESM F - erector spinae muscle force 
predicted by the model, F ML ** erector spinae muscle force derived from the RMS-EMG of the left erector 
spinae muscle, F MR = erector spinae muscle force derived from the RMS-EMG of the right erector spinae 
muscle, RAM F =* rectus abdominus muscle force predicted by the model, F RL = rectus abdominus muscle 
force derived from left rectus abdominus muscle, and F RR - rectus abdominus muscle force derived from 
right rectus abdominus muscle. Subject and task factors are explained in the text. Horizontal GRFs all 
depended on the direction of exertion, with differences in the subject factors involved. Vertical GRFs were 
affected by more subject factors. The trunk forces alt depended on the individual subject factor; the only task 
factors of importance were the direction of exertion and the cart resistance. Notice that the variabiUty in the 
measured variables was more readily explained by the selected subject and task factors than the variability in 
the predicted variables. 


tween derived and predicted means was that direction 
of exertion was not a significant factor in the model 
predicted values. Model predicted mean ESMFs were 
less than their derived counterparts, most notably for 
the subjects in the average weight category. The model 
predicted that mean RAMFs exceeded the derived 
values for light and average weight categories, but fell 
between derived values for the heavy subjects. 

Analysis of the differences between measured and model 
predicted variables 

The residuals between measured and predicted 
horizontal GRFs were affected by subject, walking 
speed, and direction of resistance (SN, WS, and DR, as 


shown in Table 5). The positive residual means for 
both horizontal and vertical GRFs reaffirm the bio- 
mechanical model’s bias towards overprediction. The 
vertical GRF residuals were affected by gender and 
height (SX and HT). The horizontal to vertical GRF 
ratio residuals were not explained by the linear statist- 
ical model. Average residual values were partitioned 
by subject body weight category for both single and 
total right leg support for the horizontal and the 
vertical GRFs. The horizontal GRF residuals were 
closest to zero for the subjects in the average weight 
category ( WT = 2) when pulling. During pushing the 
residuals stayed positive across body weight. The 
vertical GRF residuals were minimized with the light 
subjects (WT= 1) for both pushing and pulling. 
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Table 4. Average 
Weight 


trunk muscle force ANOVA summa ry partitioned by subject body weight category 
Variable K-squar e Mean Subject factors Task factors 

Trunk extensors 


Single support 

Light 

ESMF 

Light 

FMR 

Light 

FML 

Average 

ESMF 

Average 

F MR 

Average 

FML 

Heavy 

ESMF 

Heavy 

FMR 

Heavy 

FML 

Total support 

Light 

ESMF 

Light 

FMR 

Light 

FML 

Average 

ESMF 

Average 

F MR 

Average 

FML 

Heavy 

ESMF 

Heavy 

FMR 

Heavy 

FML 


0.37 

0.86 

0.95 

0.4 

0.79 

0.65 

0.38 

0.69 

0.77 


0.39 

0.87 

0.95 

0.4 

0.79 

0.67 

0.39 

0.67 

0.78 


504.5 

SN 

— 

760.4 

SN 

DR 

620.4 

SN 

DR 

534.2 

SN 

— 

1141.2 

SN 

DR 

1038.9 

SN 

DR 

622 

SN 

— 

810.9 

SN 

DR 

761.8 

SN 

WS, DR 

476.6 

SN 

OR 

748.2 

SN 

DR 

620 

SN 

DR 

481.6 

SN 

— 

1151 

SN 

DR 

1011.2 

SN 

DR 

561.3 

SN 

— 

821.8 

SN 

DR 

722.2 

SN 

WS, DR 


Trunk flexors 


Single support 

Light RAMF 0.14 

Light FRR 0-52 

Light FRL 0.4 

Average RAM F 0.24 

Average F RR 0.53 

Average FRL 0.52 

Heavy RAM F 0.32 

Heavy F RR 0.55 

Heavy FRL 0.37 


35.8 

— 

— 

3.8 

SN 

DR 

5.3 

SN 

DR 

37 

— 

— 

19.7 

SN 

DR 

7 

SN 

DR 

116 

SN 

— 

191.1 

SN 

DR 

73.8 

SN 

— 


Total support 
Light 

RAMF 

0.13 

Light 

FRR 

0.52 

Light 

FRL 

0.42 

Average 

RAMF 

0.31 

Average 

FRR 

0.51 

Average 

FRL 

0.5 

Heavy 

RAMF 

0.28 

Heavy 

FRR 

0.55 

Heavy 

FRL 

0.37 


39.1 

— 

— 

4 

SN 

DR 

5.9 

SN 

DR 

37.7 

SN 

— 

17.8 

SN 

DR 

6.5 

SN 

DR 

146.8 

SN 

— 

190.6 

SN 

DR 

67.6 

SN 

— 


Note that the predicted trunk muscle forces (ESM F and RAM F) were not dependent on the direction of 
exertion and that the model underpredicted the ESMF for the subjects of average weight. 


Residuals were formed between model predicted 
trunk muscle forces and derived left and right trunk 
muscle forces for the total right leg support phase only 
(Table 5). All of the trunk force residuals were affected 
significantly by the individual subject factor (SN). 
Residuals between predicted extensor muscle force 
and derived muscle forces for each side were large and 
negative. However, the residuals for the trunk flexors 
were smaller. Right side ESMF residuals were closest 
to zero for the heavy subjects. However, on the left side 
the ESMF residuals were closest to zero for the 
subjects in the average weight category. Right and left 
side residual averages for RAMF resembled each 


other more closely, with minimum residuals for sub- 
jects in the average weight category performing pulls. 

DISCUSSION 

The novelty of this particular model is the capability 
to calculate dynamic parameters that relate to either 
back overexertion risk or foot slip risk. Previous 
investigations of pushing and pulling tasks used static 
analysis, and predictions of internal back forces or 
ground reaction forces were not attempted. Limita- 
tions of the present model will be discussed first in 
light of the major assumptions required. The effects of 
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Table 5. Average GRF and trunk muscle force residual ANOVA summaries 


Variable 

R-square 

Mean 

Subject factors 

Task factors 

GRF residuals 

Single support 
RF y 

0.38 

12.74 

SN 

WS, DR 

RF X 

0.48 

29.19 

SX, HT 

DR 

COF 

0.09 

-0.34 

— 

— 

Total support 
RF, 

0.39 

13.15 

SN 

WS, DR 

RF, 

0.6 

31.62 

SN, SX, HT 

— 

COF 

0.09 

-0.41 

— 

— 



Trunk force 

residuals 



Total support 


ESMF( L) 

0.29 

-2785.5 

SN 

— 

ESMF( R) 

0.87 

-857.8 

SN 

DR 

RAM F{L) 

0.76 

-13.1 

SN 

DR 

RAM F(R) 

0.75 

-20.2 

SN 

DR 


The positive means for the GRF residuals indicate the model's tendency to overpredict, 
although these average residuals are small. Note that residuals were larger for the trunk 
forces. 


these assumptions on model performance will then be 
described. Finally, the validation approach used here 
will be compared with other validations of biomech- 
anical models. 

Model performance can not exceed the quality of 
the input to the model; hence the input to the model 
must be considered as a source of discrepancy between 
measured and predicted values. Sagittal plane joint 
center position measurement with the SELSPOT real- 
time position measurement system had two inherent 
problems: firstly, the location of joint centers was only 
grossly estimable with markers placed on the skin or 
clothing at an approximated joint center, as deter- 
mined by palpation. This problem has been docu- 
mented extensively before (see Zahedi et ai 1987 for 
relevant arguments) and will not be probed here. 
Secondly, the use of a lateral photodetector caused 
distortion problems due to the lens, the detector, and 
signal reflection (Gustafsson and Lanshammar, 1977). 
The steps described in the data acquisition section 
compensated for these distortions. 

Kinematic considerations played a large role in the 
determination of single vs total right leg support. 
Model GRF prediction differed depending on the type 
of support used; without the use of foot switches or a 
second force plate, the determination of left foot 
contact with the ground relied on a foot LED floor 
clearance which had a 1.0 cm tolerance and a Y 
coordinate velocity of < 10% of maximum Y foot 
velocity. Confounding this determination was the 
subject's tendency to drag the feet along the floor, 
perhaps to maximize proprioceptive cues, particularly 
when going backwards during a pull. Uncertainty in 
selecting double vs single support could cause dis- 


continuities in GRF predictions in the transition 
region. 

Several anthropometric variables were based on 
previous work done on limited samples of cadavers 
[link lengths and masses from Dempster (1955) and 
Clauser et ai (1969); radii of gyration from Chandler et 
al. (1975); diaphragm area from Morris et ai (1961) 
and Fisher (1967)]. The individual subjects performing 
the current experiments may not be described accur- 
ately by these proportionality constants and average 
areas, particularly because of the wide range of an- 
thropometries selected. Detailed individual anthropo- 
metric measures would have been required to circum- 
vent these innaccuracies, but these were beyond the 
scope of this study. 

The GRF calculations during double support came 
from moment equilibrium equations. The selection of 
the reference point for moments taken around the 
foot-ground interface was an estimate because the 
center of foot force application in the sagittal plane 
could not be determined without recording moments 
from the force platform. Therefore the proper location 
of the foot contact reference point for horizontal 
moment arms could have been anywhere between the 
heel and the toe marker (the model assumed the 
reference points were the heel marker for pushing and 
the toe marker for pulling). It was possible that 
inaccuracies in the horizontal moment arm estimates 
contributed to model prediction errors during double 
support. 

The other assumption about double support was 
that the friction utilization was the same under both 
feet. Since only one force platform was available, this 
could not be directly verified. When overcoming cart 
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inertia (or exerting push or pull forces maximally 
against a stationary object) most subjects concen- 
trated their foot forces on one foot, using the other 
foot for balance (Chaffin et ai (1983); also based on 
observation of the volitional postures and foot forces 
used by subjects in the static calibration experiments). 
This concentrated the vertical forces on a smaller 
surface area, and thus increased the contact pressure 
at that particular shoe-floor interface. How these 
static exertions relate to the dynamic case remains 
unclear. 

Since the linear regression coefficients derived dur- 
ing the static calibration experiments were used to 
transform the RMS-EMGs (obtained during the dy- 
namic experiments) into derived muscle forces for 
comparison with predicted muscle forces, they played 
a role in evaluating model performance, and so are 
discussed here. Muscle activity, as described by 
EMGs, is highly variable. Processed EMGs (i.e. integ- 
rated, average, RMS) have been shown to describe 
more accurately levels of muscle force output than raw 
EMGs, particularly in isometric situations (see 
Redfem (1989) for a complete review of these argu- 
ments). It is also well known that surface EMGs are 
susceptible to signal decrements due to adipose tissue 
interposed between the sensor and the muscle tissue 
(Basmajian and DeLuca, 1985). The subjects in these 
experiments were not of uniform somatotype, and so 
the quality of surface EMG obtained from the heavy 
subjects was suspect. Another factor may have been 
the use of isometric exertions for the calibrations. 
Most previous work correlating RMS-EMGs with 
muscle force utilized discrete isometric force levels 
only (Basmajian and DeLuca, 1985), as opposed to 
our ramp of increasing force. 

Vertical and horizontal GRFs were compared for 
both single and total right leg support because of the 
sensitivity of model predictions to the transition from 
one support phase to the other. The most rigorous 
validation of the model GRF predictions came from 
the single support phase, since none of the assump- 
tions necessary to make the double support solutions 
determinate were required. 

Measured average horizontal GRF values through- 
out the total right leg support phase were less than the 
model predictions. However, inspection of the residual 
analysis results revealed that larger intra-subject dis- 
crepancies existed when the double support phase was 
included, as opposed to those during single support. 

Performing the same comparisons for the vertical 
GRFs, single support will be discussed first. The 
amount of variability explained by the regression 
analyses was similar for measured (66%) and predic- 
ted (86%) vertical GRFs. Considering the residuals 
formed by taking the difference between predicted and 
measured vertical GRFs during single support, the 
mean average residual value was 29.2 N, indicating 
that the model overpredicted more than it under- 
predicted. 

During total right leg support, the mean average 


values were again greater for predicted than for meas- 
ured vertical GRFs. The intra-subject comparison 
disclosed that the mean average force residual for the 
total right leg support phase was essentially the same 
as for the single support phase; this implied that the 
patterns of model over- or underprediction during 
double support were not consistent from subject to 
subject, and hence were counteracted when combined 
across subjects. The statistical consensus was that 
model predictions of GRFs during double support 
were not as valid as those during single support. This 
must be attributed to both the limitations of the 
moment arms used and the assumption of the equal 
friction utilization by both feet during double support. 

Since the derived ESMFs came from trans- 
formations of measured RMS-EMG values with re- 
gression coefficients obtained from the static experi- 
ments, any limitations in the regressions will be re- 
flected in the derived ESMFs. The ANO VA analysis of 
derived ESMFs indicated that anywhere from 68 to 
74% of the variability in the data was explained by the 
independent main factors. 

Considering the analysis of residuals formed by 
subtracting the derived ESMFs from the predicted 
ESMFs , the activity of the right side erector spinae 
group was predicted better than the activity of the left 
side muscle group. On an intra-subject basis, the 
predictions underestimated average ESMFs over the 
total right leg support phase. The previous discussion 
of surface EMGs taken from varied somatotypes is 
also germane, but, given all of these qualifiers, predic- 
ted ESMF performed qualitatively the same as right 
side derived ESMF . It should be stressed that the 
support phase used throughout this study occurred on 
the right foot, and therefore it was not surprising that 
right side derived £SAf Fs seemed more realistic than 
left side values; the left side muscles were involved to 
some degree in the swing phase of the left leg, which 
occurred during right leg support. 

The analysis by subject body weight category pro- 
vided better insight into the indirect validation of 
predicted muscle forces compared with muscle forces 
derived from dynamic RMS-EMGs. Perfect corres- 
pondence of derived and predicted values would yield 
residuals equal to zero. In most cases the ESMF 
residuals were closest to zero for the medium weight 
subjects. Adding the consideration that the erector 
spinae are trunk extensors, they were expected to be 
most active during pulling. Indeed, in most cases, the 
residual parameters were closer to zero for the pulling 
tasks in this investigation. 

Derived and predicted RAMFs were examined in a 
similar manner. Based on the argument that the left 
side of the muscle group may have some involvement 
as a synergist during left leg swing, only the right side 
rectus abdominus muscle group will be discussed here. 
Minimal anterior muscle activity was seen during the 
pulling tasks. As was the case with the predicted 
erector spinae forces, more subject factors (SX, HT, 
WT) were present in the ANOVA results than for the 
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derived forces in the anterior muscles. It should be 
stressed that the linear regressions derived from the 
static calibration experiments did not, in general, 
explain as much of the variability in the RMS-EMGs 
of the rectus abdominus as for the erector spinae. This 
was attributed partly to the variation in adipose tissue 
distribution in the anterior part of the trunk as 
opposed to the posterior low back region. 

The RAM F residual parameters were closest to zero 
for either the light or the average weight subjects. The 
rectus abdominus was most active during pushing; i.e. 
the residual averages during pushing were low. How- 
ever, residual maxima (all positive) during pushing 
increased as weight increased, reflecting an increasing 
model overprediction with heavy subjects. 

Lee (1982) also reported much better model per- 
formance in predicting GRFs than EMGs. He found 
that predicted vertical GRFs correlated with meas- 
ured vertical GRFs (r 2 =0.65) and predicted horizon- 
tal GRFs correlated with measured horizontal GRFs 
(r 2 = 0.56). Mean errors were small (60 and 19 N, 
respectively). This performance was similar to that 
found in the current investigation. Lee found that 
there was a subject effect in the differences between 
predicted and measured GRFs. This was again found 
in the current investigation in the residual analysis. 

The validation approach in this investigation was 
similar to Lee (1982), with the exception that he used 
the regression relationships between isometric RMS- 
EMGs and exerted trunk forces to transform pre- 
dicted dynamic trunk muscle forces to predicted dy- 
namic RMS-EMGs, which were then compared with 
the measured dynamic RMS-EMGs. Since there have 
been no other dynamic biomechanical model analyses 
of pushing and pulling, related validation techniques 
can only be found in different tasks. Static predictions 
of mean spine compression acting at L3 have been 
correlated with mean intradiscal pressure measure- 
ments (Schultz and Andersson, 1981), with a correla- 
tion coefficient of 0.91. This same group later de- 
veloped a model to predict trunk muscle forces during 
isometric weight-holding and force resistance tasks 
(Schultz et al . , 1982), which they validated by com- 
paring their calculated trunk muscle tensions with 
mean myoelectric signal levels. Correlation coeffi- 
cients ranged from 0.34 to 0.92, depending on the 
muscle group and on the function used to predict 
muscle force. More germane to the current investiga- 
tion were predictions of dynamic trunk loading; Jager 
and Luttman (1989) developed a dynamic 19-segment 
model to assess lumbar stress during load lifting. Their 
validation consisted of comparing their model calcu- 
lations with intradiscal pressure measurements taken 
from the literature; only static holds were compared. 

Direct comparisons of predicted and measured 
ground reaction forces have been performed by others 
(Pandy and Berme, 1988, 1989). These studies simu- 
lated GRFs by assuming joint moment trajectories 
and performing open-loop (single support) or closed- 


loop (double support) analyses during normal or 
pathological walking. However, cart pushing or 
pulling tasks have not been similarly examined. 

Previous biomechanical model validation ap- 
proaches have not systematically varied subject an- 
thropometry and gender. The evidence from this 
investigation suggests that this particular dynamic 
biomechanical model was valid when a wide range of 
anthropometries was studied. Model predictions, in- 
cluding both GRFs and trunk muscle forces, were 
better during single support phases than during 
double support phases. 
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